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Abstract 

We introduce and solve exactly a family of invariant 2x2 random matrices, depend- 
ing on one parameter r\, and we show that rotational invariance and real Dyson index 
(3 are not incompatible properties. The probability density for the entries contains 
a weight function and a multiple trace-trace interaction term, which corresponds to 
the representation of the Vandermonde-squared coupling on the basis of power sums. 
As a result, the effective Dyson index (3 c r of the ensemble can take any real value 
in an interval. Two weight functions (Gaussian and non-Gaussian) are explored in 
detail and the connections with /3-ensembles of Dumitriu-Edelman and the so-called 
Poisson-Wigner crossover for the level spacing are respectively highlighted. A curi- 
ous spectral twinning between ensembles of different symmetry classes is unveiled: 
as a consequence, the identification between symmetry group (orthogonal, unitary 
or symplectic) and the exponent of the Vandermonde (f3 = 1, 2, 4) is shown to be po- 
tentially deceptive. The proposed technical tool more generically allows for designing 
actual matrix models which i) are rotationally invariant; ii) have a real Dyson index 
/? e g; iii) have a pre-assigned confining potential or alternatively level-spacing pro- 
file. The analytical results have been checked through numerical simulations with 
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an excellent agreement. Eventually, we discuss possible generalizations and further 
directions of research. 

Key words: Random Matrix, Vandermonde, correlations, Poisson-Wigner 
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1 Introduction. 

Ensembles of matrices with random elements have been widely studied since 
the pioneering works of Wigner (1) and Dyson (2) on the 'threefold way'. A 
first, gross classification of random matrix (RM) models can take into account 
i) whether the size N of the matrices in the ensemble is finite or the limit 
N — > oo is taken and ii) whether the probability distribution of the entries 
remains invariant after a rotation in the matrix space. 

The requirement of rotational invariance implies that the joint probability 
density (jpd) of the eigenvalues can be written as: 



where V{x) is a confining potential (x 2 for Gaussian ensembles) and the inter- 
action term between eigenvalues is the well-known Vandermonde determinant 
raised to the power (5. The Dyson index (3 can classically take only the values 
1,2,4 according to the number of variables needed to specify a single entry 
(1 for real, 2 for complex and 4 for quaternion numbers). This (3 index in 
turn identifies the symmetry group of the ensemble (Orthogonal, Unitary and 
Symplectic respectively). 

Thanks to the works of Mehta (3) and many others, very powerful analytical 
tools are available to deal with invariant ensembles, both for finite N and as 
N — > oo, the latter limit being usually the most interesting for RM theorists. 
However, it was very soon realized that matrices with the smallest size N = 2 
can equally well provide deep insights and trigger new ideas, the most success- 
ful one being the celebrated Wigner's surmise (3) which gives an excellent ap- 
proximation for the level spacing of bigger matrices. The study of 2 x 2 random 
matrices has since been strongly developed and it remains an active area of 
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research in mathematical physics (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) 
(16)(17). 

The purpose of the present paper is to introduce and solve exactly a family 
of 2 x 2 random matrices depending on one parameter 77. This ensemble will 
have rotational invariance but a real effective Dyson index /3 e g in an interval. 
Although it is commonly assumed that the two properties: 

• rotational invariance; 

• real Dyson index. 

are essentially incompatible, since the Dyson index of an invariant ensemble is 
strictly constrained to the values 1, 2 or 4 as described above, we will show how 
to construct explicitly a counterexample in Section 2 introducing suitable cor- 
relations among the matrix entries. The motivation for this study stems from 
two apparently unrelated issues, namely the Dumitriu-Edelman /3-ensembles 
(18) and the so-called Poisson-Wigner crossover for the level spacing (19). In 
order to make the paper self-contained, we give a brief introduction to both 
of them highlighting also the two main tasks we tackle in this paper. In sub- 
section 1.3, we provide the plan of the article. 

1.1 [3-ensembles of Dumitriu-Edelman. 

Consider the jpd (1). Does there exist a non-trivial matrix model having (1) 
as its jpd of eigenvalues for any (3 > 0? Very recently, Dumitriu and Edelman 
were able to answer this question affirmatively (18). They introduced two 
ensembles of tridiagonal N x N matrices with independent entries, whose jpd 
of eigenvalues is exactly given by (1) for general (3 > (18). These ensembles 
have been called /3-Hermite and /3-Laguerre, according to the classical weight 
their jpd contains. This result is essential for an efficient numerical sampling 
of random matrices (20) and has triggered a significant amount of further 
research (21) (22) (23) (24) (25). 

Note that the /3-ensembles, having independent non-Gaussian entries are obvi- 
ously non-invariant. Thus, the first novel task we tackle in this paper (Section 
3) is the following: 

Task 1 Design and solve exactly a (2 x 2) ensemble with: 

• rotational invariance; 

• running /3 c g > 1 ; 

• assigned classical potential (in particular, Gaussian-Hermite). 



Comments on the case /3 e fr = are given in Section 3. 
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In fact, an invariant matrix model displaying a running Dyson index would 
be of great interest: tuning the strength of the correlations between the eigen- 
values in (1) has significant importance for systems which, although endowed 
with an intrinsic invar iance, are subjected to a weak non-invariant perturba- 
tion (see e.g. (26)) and may also have important implications for lattice gas 
theory (27). Furthermore, it is a long-standing observation that nuclear sys- 
tems with two-body interactions display an average density of states whose 
profile is much closer to a Gaussian distribution (28) (29) than to a semicir- 
cle. Hence, a RM approach with the appropriate symmetries clearly requires 
much weaker, and possibly suppressed altogether, correlations among the en- 
ergy levels than those arising from (1) with integer and fixed (3. In this respect, 
the limit (3 e $ — > of our model is particularly appealing (see Section 3). 

1.2 Poisson-Wigner crossover. 

Another interesting transitional regime in quantum chaos theory, namely the 
so-called Poisson-Wigner crossover for level spacings, has attracted much at- 
tention in the past twenty years (19). In terms of the dimensionless nearest- 
neighbor spacing s, the Poisson and Wigner distributions are given by: 



and correspond to the limiting cases of classical dynamics, namely purely reg- 
ular and completely chaotic. Intermediate regimes between those two extremes 
have been intensely investigated (see (30) for a review), and interpolating phe- 
nomenological formulas have been proposed, the most famous being the Brody 
(31) and Berry-Robnik (32) distributions. The quest for a deeper understand- 
ing of such a crossover has motivated many proposals of parametrical random 
matrix models whose level spacing distribution interpolates between (2) and 
(3) (33) (34) (35) (4) (5) (6). Normally, the requirement of rotational invariance is 
the first to be dropped in those models. The reason is easy to understand: once 
this condition is imposed, the Vandermonde-coupling between the eigenvalues 
forces the level spacing P(s) to display a term of the form ~ s 13 {(3 = 1,2,4) 
and thus is very stiff, at least for small values of the gap s. No meaningful 
crossover could occur in such models for any standard choice of the confining 
potential. This problem would be overcome by an invariant model with a tun- 
able index (3 C $ and thus leads to our second unconventional task (Section 4): 

Task 2 Design and solve exactly a (2 x 2) ensemble with: 
• rotational invariance; 




(2) 



(3) 
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• running (3 e s > 0; 

• assigned level-spacing profile. 



1.3 Plan of the paper. 

The 2x2 ensemble we are going to introduce in Section 2 is completely defined 
when one assigns: 

• A symmetry group (SG) (Orthogonal, Unitary or Symplectic), correspond- 
ing to real symmetric, hermitian or quaternion self-dual matrices; 

• A weight function, to be defined in Section 2; 

• A range for the free parameter 77. 

As far as the SG is concerned, in the present study we will confine ourselves 
to hermitian (unitary invariant) matrices, although generalizations to other 
SG may be easily derived (see Section 5). The Dyson index for this Unitary 
ensemble turns out to be f3 C R = (3 — 2rj = 2 — 2t] and for this reason we call 
the ensemble r\- UE (//-Unitary Ensemble). 

In Sections 3 and 4 we make two different choices for the combination (weight 
function + range for 77) in order to tackle the tasks 1 and 2 described above. 
More precisely: 

• Section 3 : choosing as an example a standard Gaussian potential, we design 
a 77- UE ensemble which is essentially a (2 x 2) /5-Hermite model (18) plus 
rotational invariance for rj e [0, 1]. We compute analytically the marginal 
distributions of the correlated entries in subsection 3.1 and we derive ex- 
plicitly the spectral properties in 3.2. These results are then checked by 
numerical diagonalization of actual ?7-UE samples in subsection 3.3. 

• Section 4 : choosing as limiting cases the Wigner and the Poisson level- 
spacing profile, we design a 77- UE ensemble whose level spacing interpolates 
between the two cases for the parameter 77 e [1/2, 1]. Following the same 
guidelines, it is in principle possible to extend the analysis to an arbitrary 
pre-assigned level-spacing profile P^s). 

In Section 5 we discuss generalizations of this model towards different SG, dif- 
ferent weight functions and extended ranges for 77. At that stage, we will make 
comments about some emerging features of our model that appear interesting 
to be tackled in future researches. 

In Section 6, we first provide a synthetic table with a comparison of the main 
features of all the ensembles considered in this work, and then we add some 
concluding remarks. Some technical derivations are also given in the Appen- 
dices. 
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2 Main idea and the model. 



Let Pjjl/f] = P,j(^ii, • • • ,xnn) be the joint probability density of the entries 
for a N x N random matrix ensemble, depending on the parameter rj. If the 
model is required to be rotationally invariant, as in our case, two facts must 
be taken into account: 

(1) Weyl's Lemma holds (3), so Pr?|/f] can be only a function of the traces 
of the first N powers of X. We highlight this point by writing hereafter 
starred quantities (such as P* := P^[Af]) whenever they are meant to be 
written in terms of traces of powers of X. 

(2) the jpd of eigenvalues is given by: 

P IJ (Ai,...,A JV )ocP*x ni^-^ ( 4 ) 

j<k 

where the Vandermonde term comes from integrating out the 'angular' 
variables in the diagonalization X — > OAO -1 . In (4), the (3 index can take 
only the values 1, 2 or 4 according to the SG of the ensemble (Orthogonal, 
Unitary or Symplectic respectively). 

We can specialize the properties 1 and 2 to an ensemble of 2 x 2 unitary 
invariant hermitian matrices: 



X 



( „ t+is 
X V2 

y 



(5) 



where x,y,t,s are random variables taken from a jpd P v (x,y,t, s) and the 
l/v2 factors are included for later convenience. 

In this simplified case, (4) becomes: 

P r? (A 1 ,A 2 )ocP;x |A 2 — A x | 2 (6) 

and we choose to write the ?7-UE jpd of entries as: 

In (7), the weight function W* is a non-negative, normalizable and symmetric 
function of the eigenvalues, expressed in terms of the traces (TrAf, TiX 2 ). It 
may depend or not on the parameter 7] < 3/2. 

Now, we define: 

V* = 2Tr^ 2 - (TiX) 2 (8) 
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and it is easy to prove the following identity involving the rhs of (8): 

2Tr^ 2 - (TrAf) 2 = |A 2 - Ai| 2 (9) 

Through (9), the Vandermonde-squared coupling has been represented in 
terms of traces of powers of X and the jpd of eigenvalues (6) becomes: 

P,(A 1 ,A 2 )ocV^x lAs-Ail 2 " 2 " (10) 

As r\ G X C JR., the effective Dyson index (3 e n = 2 — 27/ assumes real values in 
an interval, while the ensemble keeps its rotational invariance (unlike in the 
case of /3-ensembles of Dumitriu-Edelman). The price to pay is that the entries 
are no longer independent, but get correlated through the multiple trace-trace 
interaction term V*. 

One may ask whether the crucial identity (9) is just an algebraic accident 
holding only for 2 x 2 matrices or it has a deeper origin. In fact, (9) turns out 
to be a special case of the more general identity: 

n(Aj - A fc ) 2 = detM N [X N ] = 

j<k 



N 


Si 


S2 ■ 


• Sjv-1 


si 


S2 


S3 ■ 


sn 


S2 


S 3 


s 4 ■ 


■ Sn+1 



SN-l Sn S]y + i ■ ■ ■ S 2 (N-1) 

where Sk = TrA^ (36). Since the Vandermonde-squared coupling is a sym- 
metric polynomial in the eigenvalues, it can be represented on the basis of 
power sums (37), which are nothing but traces of higher order powers of X^. 
The representation (11) precisely encodes this change of basis, which is cur- 
rently used in the context of the fractional quantum Hall effect but usually 
overlooked in RM studies. Note that the general expansion (11) can be used 
in principle to define a N x N rj-UE model, although any analytical approach 
appears very challenging in this case. 

From (10), it is clear that different choices for the weight function and the 
range X for t] can be combined to achieve a variety of results. In particular, 
we are now ready to tackle the first task described in the Introduction. 
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3 First Task: Gaussian weight function. 



Suppose we choose the confining potential to be harmonic V(x) = x 2 . It is then 
sufficient to make the simple choices W* = exp |Tr,Y 2 ) and rj G [0, 1] to 
design a /3eff-Hermite model (/3 e fr £ [0, 2]) (18) plus unitary invariance which we 
are going to solve exactly. Before doing that, we make the following important 
remark: 

Remark 1 Unlike in the case of (3-Hermite ensemble, the value /3 e fr = can be 
actually reached in 7]-UE for rj = 1. This corresponds to independent normally 
distributed eigenvalues. 

From (7) and (8) we get: 

p* _ n (in) 

v *> [2TtX 2 — (TrX) 2 } 1 ! 1 } 

where is a normalization constant. 

The resulting jpd of eigenvalues (6) can be written as: 

P,(Ai, A 2 ) = K,e-^ A ? +A ^|A 2 - Ail 2 " 2 " (13) 

where the normalization constant is given by (18): 

K,= (^r(^))" (14) 

Note that (13) is almost equivalent to a 2 x 2 /3 c fr-Hermite jpd (apart from 
the actual /3 e g = case which is not included there), although the underlying 
matrix model is very different in the two cases. For a related model with 
complex eigenvalues, see (38). 

The range of variability for 77 is largely arbitrary (see also subsection 5.3). The 
choice of [0; 1] is motivated by a nice duality between the limiting cases 77 = 
and 77 = 1 as in the following table: 





77 = 


n = 1 


Correlation among Eigenvalues 


Strong 


Absent 


Correlation among Entries 


Absent 


Strong 



Table 1 

The jpd of entries factorizes for 77 = and gives rise to GUE (Gaussian Unitary 
Ensemble) with strongly correlated eigenvalues. On the contrary, for 77 = 1 the jpd 
of eigenvalues factorizes and the eigenvalues become i.i.d. normal variables, while 
the entries are strongly correlated. 
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Although (12) defines completely our ?7-UE model, it is instructive to com- 
pute analytically the marginal distributions for the set of correlated variables 
(x, y, t, s) for two reasons: 

(1) These results are numerically implemented in subsection 3.3 to generate 
and diagonalize actual samples of rj- UE matrices. The numerical results 
will be compared with the spectral properties derived in subsection 3.2. 

(2) The marginal distributions deviate smoothly from the GUE factorized 
marginals as r\ departs from zero, and thus they provide quantitative 
information about the onset of correlations between the entries. 



3.1 Marginal distribution of entries. 



From (5) one has: 



TiX = x + y (15) 
Tr* 2 = x 2 + y 2 + t 2 + s 2 (16) 



so that (12) implies: 



e 2 v 

P v (x, y, t, s) := C„ [2 ( x2 + y2 + t2 + s2) _ {x + y)2]ri (17) 

The first task is computing the normalization constant C v , for which the fol- 
lowing integral is needed [(39) formula 3.382(4)]: 

where T(x,y) is the incomplete Gamma function. 

The constant C" 1 is given by: 

x r°° r°° dx dy dt ds e -\( x2 +y 2 +t 2 +s 2 ) 
r > = loo " ' Loo [2(x 2 + y 2 + t 2 + s 2 ) -(x + y) 2 ]v 

which becomes upon the change to polar coordinates (t, s) — > (r, 9): 

C' 1 = 2ix r r dx dy e-^ x2+v2) X( V , x - y) (20) 

J—oo J—oo 

After further simplifications and the change of variables (x + y,x — y) — > (<f),uj) 
the double integral decouples and we get: 

C" 1 = vr 3 / 2 2 3 - 2 T (J - 77) (21) 
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Now, we can compute the marginal distribution for the variable x following 
essentially the same steps (hereafter we use the notation pfc j7? with k = 1,2,3 
to denote the marginals with k variables): 



/oo roo 
■■■ dy dt ds P v (x,y,t,s) = 
-oo J — OO 
r\ ex /"OO 

= 2nC v — J o dfij ^ dy e-^'+^Jfo, y - 



2 2 - 2 %C„j- / duTil-r]^ 2 ) d^ 

CLX J— co J uj 

/oo 
dcuT(l-r],uj 2 )e 
-oo 



(fy * 9) (x) 

2^r (I - n 



(22) 



where in the last passage we have employed the symmetry uj — > —uj of the 
integral and the convolution (*) is between the functions: 



f v {y) :=r(l -r,,y 2 ) g(y) := e"» (23) 



In the limiting case rj — > 0, we expect to recover the pure GUE marginal 
distribution for the entry x, which is simply a standard Gaussian. Taken into 
account that f v =o(y) = g(y) and T(3/2) = y/n/2, one gets: 



p -x 2 /2 

due-^e'^-^ 2 = —== (24) 



Pi,„=oW = - 

vr J-00 V27T 



as it should. 



A careful asymptotic analysis (see Appendix A) of the convolution integral 
(22) gives for x — > ±00: 



P l„(g)« r- 1 (-) 2 \~ x2/2 (25) 

* M ! 2^(3/2 -rj) \2J V ; 



and for 77 > the decay is faster than Gaussian due to the power-law prefactor, 
in agreement with the plots in figure 1. 
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Pi, 77 ( x ) 




X 



Fig. 1. Plot of the marginal distribution pi jr) (x) for different values of r\ [Orange 
0.8; Green 0.5; Blue 0.2]. The limiting Gaussian distribution r/ = is also plotted 
[Violet]. 



The other marginals can be computed as well: 



P2,r,(x,y) = 
P3, v (x, y,t) 



e 4 



r 1-77, 



[x - yf 



e -\{x 2 +y 2 +t 2 ) 
7r 3/2 2 3-2^ r (| _ ^ 



gV^yg2T(77 ~ 1/2) 1^(1/2, 3/2 - 77, a/4) 

rfo) 

+ 2 1 / 2 ~ 2 T(l/2 - 77) ^(77, 1/2 + 77, a/4)} 



(26) 



(27) 



where we defined a = (x — y) 2 + 2t 2 and i-Fi(&, 0) is the Kummer confluent 
hypergeometric function. 

At this stage, we make the following important remarks: 



(1) For — ^ 0, it is straightforward to check that (26) and (27) reproduce 
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the expected GUE factorized distributions: 




(28) 



2tt 



P3 1?7 =o(^,2/,0 



(2tt) 3 / 2 



(29) 



(2) For ?7 — > 1/2 (a noticeable special case, see next subsection), all the 
marginal distributions remain well-defined. In particular, the apparent 
divergences of the Gamma functions in (27) cancel out and the final 
density reads: 



where K (x) is a Modified Bessel Function of degree of the Second 
Kind. 



3.2 Spectral properties. 

As already discussed at the end of the previous subsection, the case T) — 1/2 
is particularly interesting, as the jpd of eigenvalues (13) collapses onto the 
Gaussian Orthogonal Ensemble (GOE) one (f3 e g = 1). The matrices belonging 
to GOE have independent and real entries and the orthogonal group as SG. 
Instead, 77-UE reproduces the GOE spectral statistics while having complex 
and correlated entries and the unitary group as SG. This is a first example of 
a curious phenomenon, which we call spectral twinning (see Section 5.1) be- 
tween ensembles having the same spectral properties (same jpd of eigenvalues) 
but different SG (different number of independent real variables). As a con- 
sequence, the connection between the exponent of the Vandermonde /3 e ff and 
the SG of the ensemble becomes much more blurred and potentially deceptive 
than for the classical f3 — 1, 2, 4 ensembles. 

In this section and in Appendix B, we compute for completeness the average 
spectral density for our rj- UE model with Gaussian weight function, as this 
calculation does not appear to have been carried out explicitly before. The 
raw level spacing has already been computed in the context of a 2 x 2 (3- 
Hermite ensemble in (40). 

The average spectral density p v (X) and the gap probability P v (s) are given as 



P3,r,=l/2(x,V,t) 



1 



e a/8 K (a/8) 



(30) 



7T 3 / 2 2 5 /2 
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usual by: 

/oo 
P r) {X,X 2 )dX 2 (31) 
-oo 
roo 

P v (s) = / P r ,{X 1 ,X 1 + s)dX l (32) 



where in (32) and hereafter, s is meant as the raw spacing, without any un- 
folding procedure being performed on the spectrum. 

They can be computed exactly from (13) as (see Appendix B): 

-A 2 

PvW = 2 3/2-„ ~ ^ V2, A 2 /2) (33) 
P v (s) = 2v^FK, s 2 - 2 "e- s2 / 4 (34) 

where is defined in (14). 

It is easy to check that (33) recovers for rj = 0, 1/2, 1 the expected spectral 
densities (GUE, GOE and purely Gaussian, respectively): 

Po(A) = l il AJ (35) 

Zy In 



e 



-A 



2 (2 + V2^e A2 / 2 A erf(A/ v / 2)) 



Pi/ 2 (A) = ^ (36) 

p -A 2 /2 

Pi (A) = -= (37) 

V Z7T 

where erf(z) = (2/^) / 2 e-' 2 dt. 

The results (33) and (34) are numerically checked by actual diagonalization 
of ?7-UE matrices in the next Section. 



3.3 Numerical Simulations. 



We report in this section the results for the spectral density and the level 
spacing, obtained by direct sampling of 77- UE matrices. 

The algorithm proceeds as follows: by rejection sampling (41) we draw a ran- 
dom number x from the marginal distribution pi j7? (x) (22). Then, from the 
marginal distribution (26), we determine the conditional probability: 

Pv(y\ x ) = J ( ^ ( 38 ) 
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and again we draw a random number y from (38). This procedure is iter- 
ated through the higher order marginals up to identifying the four variables 
(x, y, i, s) from which one sample of 77-UE is constructed. Each sample is then 
diagonalized and we give a histogram of the eigenvalues and of the gaps be- 
tween the two eigenvalues over a total number of 75000 samples for each plot. 

We include three plots for the spectral density (Fig. 2,3,4) (77 = 0, 0.45, 0.75 re- 
spectively) and three plots for the gap probability (Fig. 5,6,7 ) (77 = 0, 0.5, 0.75 
respectively). On top of the histograms, we plot the theoretical results (33) 
and (34). The agreement is excellent. 

0.35 



0.3 

0.25 
0.2 

Q. 

0.15 
0.1 
0.05 



-6 -4 -2 2 4 6 

X 

Fig. 2. Plot of the average density of states Prj(^) for r] = 



4 Second Task: Non-Gaussian weight function. 

In this section, we show how to devise a weight function depending on the 
parameter 77 e [1/2, 1] such that the level spacing for the corresponding 77- UE 
ensemble develops a Wigner-Poisson crossover. In principle, the solution we 
offer can be taken as a guideline for the general problem described in the Intro- 
duction as Task 2. Even though the two above cases (Poisson and Wigner) ab- 
sorb the vast majority of literature on these issues, nevertheless few instances 
of 'non-standard' gap distributions have been also reported (42) (43) (44). Our 
hope is that the tool we propose here may be used to comprise even these fairly 
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0.35 



0.3 



0.25 - 



0.15 



0.05 



o histogram 
— theoretical 




Fig. 3. Plot of the average density of states Pr?(A) for rj = 0.45 



0.35 



0.25 - 




0.15 



0.05 



Fig. 4. Plot of the average density of states p v (X) for rj = 0.75. Beyond the value 
r] = 0.5, the density becomes unimodal. 

anomalous cases into the universal and otherwise very successful framework 
of RM theory. 
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0.45 r 



0.35 



o histogram 
— theoretical 



0. 




0.15 - 



0.05 



Fig. 5. Plot of the gap probability P v (s) for r\ = 




s = \ 2 -X, 



Fig. 6. Plot of the gap probability Prj(s) for rj = 0.5 
The starting point is the general jpd of eigenvalues (10): 

P, (vv) (Ai, A 2 ) oc W;(X U A 2 )|A 2 - Axl 2 " 2 " (39) 

where the superscript (W) recalls that the weight function is still to be deter- 
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0.5 



o histogram 
theoretical 




Fig. 7. Plot of the gap probability -Pjj(s) for 77 = 0.75 



mined. 



The gap probability for a general weight function can be formally written as: 

/oo roo 
/ rfA 1 rfA 2 W,;(A 1 ,A 2 )|A 2 -A 1 | 2 - 2 ''x 
-00 J— 00 

/oo 
W*(Ai,Ai + s)dAi (40) 
-00 

Apart from being symmetric in the eigenvalues 2 , non-negative everywhere 
and normalizable, the sought W*(Ai, A 2 ) should satisfy the following con- 
straints: 



/oo 
W* =1 (£, x + s)dx oc exp(— \s\) 
-00 

W* =1 / 2 (a;, x + s)<£c oc exp(— as 2 ) 



(41) 
(42) 



where a is a numerical constant 3 , in such a way that (40) reduces exactly to 
Poisson for 77 = 1 or Wigner for 77 = 1/2. Indeed, the Wigner distribution in (3) 
corresponds to the gap probability for a GOE ensemble of 2 x 2 real matrices 
(Aiff = 1) and thus has to be realized in 77- UE for 77 = 1/2 and for a Gaussian 



2 This is simply because the joint distribution of eigenvalues can not depend on 
how one labels the eigenvalues. 

3 The explicit value for a is 1 if the spectrum has not been unfolded and 7r/4 in 
the other case. However, specifying a is not crucial for what follows. 
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weight function (42). The other limit is when there are no correlations among 
the eigenvalues (77 = 1) and thus one may expect a Poisson distribution for 
the level spacing (41). 

Note that the prefactors in (41) and (42) can be restored at the end by nor- 
malization and in Eq. (41) we use the absolute value of s since the function 
must be symmetric in s. 

In order to find an appropriate weight function satisfying the given constraints, 
we first make the ansatz 4 : 

w;(x,y) = Mz)<i>M ( 43 ) 

where (j) v (x) is an even function of x. 
Now, from (41) one has the convolution: 

4>i(x)(j)i(s — x)dx = exp(— \s\) (44) 

which gives in Fourier space: 

Thus: 



Inverting (46), we get: 



Uk) = ^ (46) 
01 (1) = — K (x) (47) 

7T 

where Kq(x) is a Modified Bessel Function of degree of the Second Kind. 
Similarly, for rj = 1/2 one gets: 

<P l/2 {x) = [ia/n] 1 ^- 2 ^ (48) 

In order to obtain the function interpolating between (47) and (48), we notice 
that Kq(x) has the following integral representation (39): 

K (x) = - / dt . (49) 

2 Jo t 

where 7 is any real and positive parameter. 



4 Other non-factorized weight functions may exist. 
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Then one can write: 

r°° dt ( x 2 \ . . 

««) = BlJ £ T exp(- 7( --j (50) 

where B x = l/(izy2). On the other hand, 01/2(2?) can also be written trivially 
in a similar integral representation as: 



>l/2 



= Bi/2 exp (-7^ - x 2 /4^ (51) 



where B 1/2 = 2- 5/2 a- 3/4 vr- 1/4 and 7 := 1 /(8a) . Thus, for general 1/2 < 77 < 1, 
a natural definition of ^(x) interpolating the two bordering cases would be: 



(52) 



The corresponding weight function W*(x, ?/) is given by the product in (43) 
and satisfies all the given constraints. We will refer to this weight function as 
a Generalized Bessel weight. 

The reader may be puzzled by the non-standard expression (52) and may 
wonder whether the resulting weight function W*(Ai, A2) is indeed expressible 
in terms of traces of powers of X, a fact which is not obvious at first sight. 
Actually, this can be shown easily by expressing the individual eigenvalues as: 



si ± \/2s 2 — s? 
Ai, 2 = V (53) 



where Sj = TrX^. 



The gap probability can then be computed from (40). For arbitrary rj, the 
integral in (40) is difficult to carry out for all s. However, the large s tail of 
P( w '(s) can be easily derived by the saddle point method and has the following 
behavior: 

Pi m (s) oc exp [—77(277 - 1)- 1+1 / 2 " (2a)- 1+1/v s 1 ^] . (54) 

One can check easily that it reduces to the known Wigner and Poisson cases 
for rj = 1/2 and 77 = 1 respectively. Thus our ensemble smoothly interpolates 
the tail of the gap distribution between the two known cases, although the full 
behavior of P^ w \s) is different from any previous proposals (e.g. Brody and 
Berry- Robnik distributions). 

Also the spectral density can be numerically investigated after some interme- 
diate algebraic steps that we include in Appendix C for completeness. The 
resulting plots for different values of ( := 277 — 1 (0 < ( < 1) and a = 1 are 
reported in Fig. 8. 
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Fig. 8. Plot of the spectral density p^(X) for the following different values of £: 
(dotted green), 0.45 (long-dashed blue), 0.9 (short-dashed magenta). Note the pure 
GOE density for £ = and the peculiar trend towards an integrable divergence at 
the origin as £ — > l - . 

In summary, the Generalized Bessel weight function W*(Ai, A 2 ) = r) (Ai)0 r) (A 2 ), 
where (f> n {x) is given by (52) and Ax,2 are expressed in terms of traces as (53) 
generates a 77- UE model having: 

• rotational invariance; 

• < [3 eS < 1; 

• prescribed level spacing (interpolating between Poisson and Wigner); 

• a novel transitional profile for the spectral density as documented in Fig. 8. 

and is to be regarded as complementary to the model proposed in (6). Note also 
that the above results are still valid in the range [0, 1/2] for 77, corresponding 
to a GUE-Poisson crossover. 



5 Generalizations. 

The work presented here can be extended in several directions. We would like 
to offer a list of issues that can be tackled in future researches. 
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5.1 Different Dyson class (3. 



We confined our investigation to hermitian matrices {(3 = 2), since the identity 
(9) involves exactly the exponent 2 for the Vandermonde, but it is not harmful 
to consider the following obvious relations instead: 



|A 2 - Ai| = ^2TyX 2 - (TrX)' 



I A2 — Ai I 



2Tr^ 2 - {TrX) 



(55) 
(56) 



and reformulate the model for real symmetric and quaternion self-dual Her- 
mitian matrices respectively: 



X-u 




Xm 





X 





t + is 


u + iv 







X 


—u + iv 


t — is 


t 


— is 


—u — iv 


y 





u 


— iv 


t + is 





y 



(57) 



having respectively 3 and 6 real independent variables. 

Example: real symmetric matrices with Gaussian weight function. In this 
case, the jpd of eigenvalues reads: 



1 



Pf)(\i, A 2 ) cx exp --{Xf + X>) |A 2 - Ai 



(58 



where we rename the running parameter as < fj < 1, and the jpd of entries 
is given by: 



lT r v2 

e 2 



P ^ = C « [2TrXl - (TrX R )*]W (59) 
The expression (58), when compared with the corresponding one for hermitian 
matrices (13), leads immediately to the following observation: if the parame- 
ter rj is chosen in [1/2, 1] and 17 := 2t] — 1, then the two ensembles (real and 
complex) get twinned, i.e. they share the same spectral properties, despite 
belonging to different classes of invariance and having even a different num- 
ber of independent variables. This spectral twinning is a curious byproduct 
of our construction, which was already remarked in subsection 3.2. While it 
is premature to imagine possible physical application for this, nonetheless we 
believe that this peculiar property, which does not hold for any classical in- 
variant ensemble and of course neither for the /3-ensembles, deserves further 
investigations and may be related to some group-theoretical features of 77-XE 
(X=0,U,S) yet undiscovered. Note also that this property would survive even 
for N x N 77-XE and holds for any acceptable weight function. 
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5.2 Different classical weights. 

A whole group of novel 2x2 ensembles can be generated by choosing different 
weight functions among the classical ones, as: 

(1) Laguerre and Jacobi (to make contact with (18)); 

(2) fixed and restricted trace ensembles (45) (46) (47); 

(3) quartic and higher order potentials (48); 

(4) power laws (49) (50) (51); 

In view of subsection 5.1, all the above can be generated starting from real, 
complex and quaternion entries and a number of twinnings can be found. In 
particular, the /3-version of (3) and (4) has not been constructed as an actual 
random matrix ensemble so far, while this problem can be tackled using the 
method we presented here. 



5.3 Extended range for r\. 

Provided that the weight function decays fast enough at infinity in order to 
ensure convergence of the integrals involved, it is possible to extend consider- 
ably the range of variability for rj in the two examples above as well as in any 
future study. Taken the Gaussian case as an example (Section 3), the following 
extensions can be considered: 

• r\ < 0: a negative rj in (13) enhances (instead of suppressing) the correlations 
among eigenvalues and extends the range of variability for the effective 
Dyson index from [0, 2] to [0, +oo). 

• l<?7<3/2: this case is even more interesting as it allows to generate an 
Anti-?7-UE ensemble whose jpd of eigenvalues would be given by: 



i.e. with a negative Dyson index. The most immediate consequence is that 
the peculiar level repulsion has to be replaced by a fairly uncommon (at 
least in RM studies) level attraction. The tendency of energy levels to cluster 
instead of repelling apart has been found in several disordered many-body 
systems (52) (53) (54) (55) but has not received equal attention in the context 
of invariant RM. The Anti-77-UE is likely to lead to a 'non-Wigner'-surmise 
for the level spacing to be compared with the studies cited above. We offer 
this idea as our last contribution in this paper. 
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6 Conclusions. 



Before summarizing the main results of this paper, we propose the synthetic 
Table II, containing the most relevant features of the ensembles considered in 
this work. 



Name 


Weight 


Inv. 


Indep. 


Entries 


Size 




/3-Hermite 


Gaussian 


N 


Y 


R 


N > 2 


> 


GUE 


Gaussian 


Y 


Y 


C 


N > 2 


2 


GOE 


Gaussian 


Y 


Y 


R 


N > 2 


1 


77-UE 


Gaussian 


Y 


N 


C 


N = 2 


[0,2] 


Generalized Bessel 


Y 


N 


C 


N = 2 


[0,1] 


Anti-77-UE 


Gaussian 


Y 


N 


c 


N = 2 


[-1,0] 


/3-XE 5 


Arbitrary 6 


Y 


N 


R,C,H 


N = 2 


Arbitrary 7 



Table 2 



Columnwise: Name of the ensemble, Classical weight or weight function, Rotational 
Invariance, Independent Entries, Type of Entries, Size of the matrices, Range for (3 e g 
(exponent of the Vandermonde) . In the last row, different combinations of symmetry 
class, weight function and range for r\ may be exploited according to Section 5. 

The main results of the paper can be summarized as follows: 

(1) Although the /3-index of an invariant ensemble is determined uniquely by 

its symmetry group, through a rather unusual expansion of the Vandermonde- 
squared on the basis of power sums it is possible to neutralize (entirely 
or partially) the coupling between the eigenvalues introducing suitable 
correlations among the entries. Hence, one has to be careful in deducing 
the SG of the ensemble from the exponent of the Vandermonde coupling 
/3 e ff, as this connection may be deceptive. 

(2) Using this tool, we have constructed an invariant 2x2 version of the 
/3-ensembles of Dumitriu and Edelman. This matrix model is completely 
defined assigning a Symmetry Group (Orthogonal, Unitary, Symplectic), 
a weight function and a certain range for 77: 

• For the Unitary case, with Gaussian weight function and < rj < 1, 
both the marginal distribution of the entries and the spectral proper- 
ties have been computed analytically and have been tested by numerical 

3 See Section 5.1 

4 See Section 5.2 

5 See Section 5.3 
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sampling of 77-UE matrices. The case (3 c g = 0, corresponding to inde- 
pendent normally distributed eigenvalues, is particularly interesting and 
can be obtained for 77 = 1. 
• For the Unitary case, with a Generalized Bessel weight function and 
1/2 < 77 < 1, we generate 77-UE matrices with a level spacing profile 
interpolating between Poisson and Wigner. Both analytical and numer- 
ical results have been provided. 

(3) Unlike the classical invariant ensembles, our ensembles belonging to dif- 
ferent symmetry classes may display the same spectral properties for a 
suitable range of the free parameter. We called this curious phenomenon 
spectral twinning and we leave a deeper understanding of it as an open 
problem. 

(4) An extended range for 77 may lead to invariant ensembles whose eigen- 
values tend to cluster instead of repelling apart due to a negative Dyson 
index (3. We are not aware of any previous proposal in this direction, 
even though this behavior is fairly common in the study of disordered 
many-body systems. 

While generalizations to bigger sizes N > 2 appear difficult to be tackled 
analytically, nonetheless the model we presented here displays non-trivial and 
often surprising features, which make us hope that the proposed technical tool 
may prove useful in future RM studies. 
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A Asymptotic Analysis of p l r? (x). 

We start from the exact expression of the marginal distribution (22): 
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We first divide the integral into 2 parts: over [— oo,0] and [0, oo]. In the first 
part, we also make a change of variable y — > —y. This gives: 



Pl '" W 2^(3/2-7?) 

+ / dyT(l- V ,y 2 )e-^ 



/ dyT(l- V ,y 2 )e-^ + 
Jo 



2v/?r(3/2 - 77) 



[Ji(x) + J 2 (x)] 



(A.2) 



Consider first J±(x). Clearly, the most important contribution to the integral 
comes from the region around y = x for large x. Since y is large in this regime, 
we can replace T(l — r],y 2 ) by its leading asymptotic form: 

r(l-r,,y 2 )*e-y 2 y- 2 *. (A.3) 
Substituting this in the expression for Ji(x) we get, to leading order: 

Ji(x)w !™ dye- y2 y-^e'^ 2 (A.4) 

J a 

where the lower limit a ~ 0(1). Rewriting the exponential we get: 

J x (x) « e~ x2 l 2 r dyy- 2r >e- 2 (y~ x W 2 (A.5) 

J a 

Making a change of variable, \[2 (y — x/2) = u, we get for large x: 

J^x) ^ ^e- x /2 {x/2)- 2r ' due~ u (A.6) 



y/2 JV2(a-x/2) 

In the x —>■ 00 limit, the lower limit of the integral can be replaced by —00. 
This finally gives: 

J x (x)K^e-*l*{x/2)-*> (A.7) 

Now for the integral J2(x), for large x, the most important contribution comes 
from the region y — 0. Thus, it is easy to see that for large x, J%{x) ~ 
~ Vi V ) — x where b v is a constant. Thus, clearly J 2 (x) 
decays faster than J\{x) and hence can be neglected for large x. 

Thus, for large x, to leading order we get: 



1 /x\' 2r ' 



e~ x2/2 (A.8) 



2 v / 2r(3/2-r/) V2, 
It is easy to check that for 77 = 0, (A.8) reduces to the standard Gaussian. 
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B Spectral Density for the Gaussian weight function. 

Consider the joint distribution of the two eigenvalues (13): 

P,(A 1? A 2 ) = Kie-M+W 2 |A 2 - Ail 2 - 2 " (B.l) 
where the normalization constant K v is given by (14). 

We want to compute the average density of states which is precisely the 
marginal distribution: 

/oo 
P„(A,A 2 )dA 2 (B.2) 
-oo 

Substituting the jpd from (B.l) into (B.2) and making a change of variable 
A 2 — A = x we get: 

Pv (\) = K„ e" A2 / 2 [°° e~^ 2 / 2 \x\ 2 ~^ dx (B.3) 

J — oo 

Divide the integral into two parts: over [—00, 0] and [0, 00] and write it as: 
p v (X) = K v e- x2 ' 2 f r [e^ 2 ' 2 + e"^ 2 / 2 } x 

.J 

x x 2 - 2r > dx] = Kr, e" A2/2 [Ji(A) + J 2 (A)] (B.4) 
The first integral can be rewritten as: 

POO 

IJX) = e~ x2 / e^ 2 / 2 x 2 - 2 ^x = e" A2 / 4 r(3 - 2 V ) x 
Jo 

x J D 2 ,_ 3 (-A) (B.5) 

where D p (z) is the parabolic cylinder function of index p and argument z (39). 
Similarly I 2 (A) = A). Adding the two we get: 

p„(X) = K, T(3 - 277) e~ 3A2 / 4 [D 2v _ 3 (X) + D 2 ^ 3 (-X)} (B.6) 
We next use the identity (9.240) of (39): 



2^2+1^-^/4 

r((l-p)/2) 



D p (z) + D p (-z) = — -y. ° l F 1 (-p/2, 1/2, ^ 2 /2) (B.7) 



where iFi(a, (3, z) is the Kummer confluent hypergeometric function defined 
by the series: 



a z a(a + 1) z 2 



Ma^tz) = 1 + -— + ~\ - -r + ... (B.8) 
U ; + 1) 2! 1 ' 
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Using this identity in (B.6) we get: 



K,2^/ 2 0Fr(3-2r,) , 
P*W- T{2 _ v) e x 

x ^(3/2-77,1/2^72) (B.9) 



We then use the explicit expression of from (14) to get: 

"' (A) = rW2-Sp-l) e ' A ' lFl(3/2 -''- 1/2 - A ' /2) (R10) 

We can further simplify the constant term by using the Gamma function 
identity (doubling formula): 

T{2z) = ^^T{z)T{z + 1/2) (B.ll) 
V 71 " 

This then gives our final formula: 



-A 



2 



PvW = 23/2-,^ i F i( 3 / 2 - ^ V2, A72) (B.12) 



C Spectral Density for the Generalized Bessel weight function. 

The spectral density for the Generalized Bessel weight function (43) is given 
by: 

/oo 
dyftMIA-i/l 1 -' (C.l) 
-00 

where the proportionality constant is just given by the normalization condition 
and we have slightly changed the notation <p v (x) — > (j>^(x), where: 

(j)J x ) := [°° dt t~ c exp ( — 2axH' c ) (C.2) 

Jo V 8a / 

and < C < 1. 

We aim to rewrite (C.l) in a form suitable for a numerical implementation, so 
from now on we will drop all the constant in front and we will normalize the 
numerical density to 1 at the very end. 

First, we define: 

;>oo 

FAx):= dss^Ms + x) (C.3) 
Jo 

and it is easy to check that: 

P? (A)oc0 c (A)[F c (A) + F ? (-A)] (C.4) 
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Exploiting (C.2) and changing the order of integration, F^(x) can be rewritten 
as: 

F{{x) oc / dt t~^ exp I 2ax 2 t~ < * I / ds s 1 ^^ exp ( — ws 2 — ns) 

8a ° (C.5) 

where w := 2at ^ and k := 4axi 



The s-integral can be carried out analytically and thanks again to the identity 
(9.240) of (39) we get: 

FAX) + F c (-X) oc /°°dt r^exp (-— - 2aX 2 t~ c ) x 
jo V 8a / 

x x Fi (l -C/2, l/2,2aA 2 t- c ) (C.6) 



Hence, the spectral density (C.4) is determined as the product of two integrals 
that are easily evaluated numerically. The results for different values of ( and 
a = 1 are plotted in Fig. 8. 
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